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Abstract 



We describe a mechanism of the primordial black holes formation that can explain the existence of a 
population of supermassive black holes in galactic bulges. The mechanism is based on the formation of 
black holes from closed domain walls. The origin of such domain walls could be a result of the evolution 
of an effectively massless scalar field during inflation. The initial non-equilibrium distribution of the 
scalar field imposed by background de-Sitter fluctuations gives rise to the spectrum of black holes, 
which covers a wide range of masses - from superheavy ones down to deeply subsolar. The primordial 
black holes of smaller masses are concentrated around the most massive ones within a fractal-like 
cluster. 
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1 Introduction 



There is growing evidence for tiie existence of massive dark objects at tlie centers of galaxies, 
which are presumably supermassive black holes (SMBHs) of mass IO^-IO^^Mq. The existence of 
SMBHs was predicted by arguments based on the active galactic nuclei (AGNs) and quasistcllar 
objects (QSOs) energetics and demography [1, 2, 3]. Also a variety of techniques exists in the 
mean time to detect and weigh SMBHs, using measurements of stellar or gaseous kinematics [3]. 
This has revealed that many, and may be all, quiescent galaxies host SMBHs as well [4], 
including the Milky Way [5]. The mass of a SMBH is always of the order of 0.1% of the galaxy 
bulge/spheroid mass [3]. An even better correlation exists between the mass of SMBHs and 
the stellar velocity dispersion in the host galaxy, Mbh oc ct^ [6], which holds 3-4 orders of 
magnitude. 

The origin of SMBHs is far from certain, but several theories have been advanced. The 
above-mentioned close relationship between BH masses and the stellar velocities dispersion 
may be interpreted in various ways, and ultimately one is confronted with the problem of what 
came first, the host galaxy or the SMBHs. It is not inconceivable that SMBHs are purely the 
result of the galactic dynamics and their merging history; yet no one has proposed a concrete 
mechanism for converting stellar mass objects into others 6-10 orders of magnitude larger. 
In particular the paradigm of the formation of SMBHs as a pure result of galactic dynamics 
is hard to match with the discovery of quasars at redshift > 6 [7]. The point is that the 
formation of SMBHs through the growth of BH seeds by gas accretion is supported by the 
consistency between the total density in quasars light and the SMBH mass density in local 
galaxies, adopting a reasonable accretion rest-mass-to-energy conversion efficiency [8]. From 
the other side, SMBHs must be present by z > 6 to power the above-mentioned QSOs. It has 
been argued [9] that if they grew by accretion from smaller seeds, substantial seeds of mass 
> IO^Mq must have been present already at ^ ~ 9 to have sufficient time to grow up to a 
typical quasar black hole mass ^ 1O^M0. Thus a very massive progenitor at 2; ^ 9 is very 
likely. At that epoch large galaxies do not yet exist. In the hierarchical formation scenario, 
they assemble via infall and merging of small clumps of baryonic gas and dark matter. Under 
such circumstances it seems that in order to form a huge BH in a gasodynamical process, two 
generic conditions must be satisfied, namely the gas must be strongly self-gravitating; and 
star formation must be inefficient, allowing a large supply of cold gas to be converted to a 
BH. The second constraint is of particular importance, since star formation is a competitor 
of black hole formation. So it is clear that to prevent the formation of stars, small-scale 
density perturbations should be effectively washed out or damped. A way to achieve this is not 
trivial. Indeed the power CDM spectrum for lightest supersymmetric particles, for instance, 
at subgalactic scales [10] implies that at 2; > 9 all fluctuations below the mass scale IO^Mq 
went to the non-linear regime, requiring some special mechanisms to dilute small scale density 
perturbations and consequently suppress star formation processes [9, 11]. So, summing up the 
above arguments it seems appealing that central primordial objects of mass > IO^Mq preceded 
the formation of the first galactic bulges and the beginning of quasar activity. 

Previously [12, 13, 14] we proposed a mechanism for the formation of the massive primor- 
dial BHs (PBHs), which opens a possibility that primordial objects of the mass > 10^ Mq 
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were created before any galaxy formation processes started to develop. The mechanism is 
based on the fact that BHs can be created as a result of a collapse of closed vacuum walls 
formed during a second-order phase transition with strongly non-equilibrium initial conditions. 
The required initial conditions get generated in inflationary dynamics of a scalar field, with 
a potential possessing an effectively flat valley [12]. A wide class of models possessing that 
symmetry-breaking pattern can be effectively described by a pseudo-Nambu-Goldstone (PNG) 
field and corresponds to the formation of an unstable topological defect in the early Universe. In 
particular, PNG fields received substantial attention in connection with curvaton proposal [15]. 

The mechanism [12, 13, 14] assumes the succession of phase transitions, related with the 
consequence of spontaneous and manifest U(l) symmetry breaking. The crucial point is that in 
the considered case spontaneous symmetry breaking takes place during inflation so that neither 
angular nor radial fluctuations can lead to formation of strings. It makes closed walls the only 
topological defects, which are formed after manifest breaking of the underlying U(l) symmetry. 
For a wide range of realistic parameters of the considered model closed walls turn out to be 
stable relative to the formation of holes bounded by strings in them. The distribution of walls 
is also practically not changed by wall collisions. It makes possible to deduce the spectrum and 
space distribution of BHs from the distribution of closed walls formed after the manifest U(l) 
symmetry breaking. 

Actually it has long been known [16, 17, 18, 19, 20, 21, 22, 23] that the PBH formation 
is possible in the early Universe. PBH can form if the cosmological expansion of a mass 
stops within its gravitational radius, which happens when the density fluctuation exceeds unity 
on a scale between the Jeans length and the horizon size. Other possibilities are related to 
the dynamics of various topological defects such as the collapse of cosmic strings [19] from 
the thermal second-order phase transition or to the collisions of the bubble walls [20, 21] 
created at the flrst-order phase transitions. Formally, there is no hmit on the mass of a PBH 
that forms when a highly overdense region collapses: it is only needed to form, at inflation, 
the spectrum of initial density fluctuations that provides the proper inhomogeneity at the 
desired scale [24]. From the other side, to preserve the successful predictions of the Big Bang 
nucleosynthesis (BBN) and general features of the standard Big Bang scenario, conflrmed by 
observations (see [23] for a review), there is no room for signiflcant cosmological inhomogeneity 
once nucleosynthesis begins; the PBHs therefore seem not to have sufficiently large masses, 
being more likely to represent the specific form of a modern collisionless dark matter, rather 
than collecting into a large BH, which could play the role of a seed for the galactic nuclei. 

The mechanism described in [12] offers a way around this, allowing the formation of closed 
vacuum walls of a size significantly exceeding the cosmological horizon at the moment of their 
formation. At some instant after crossing the horizon, such walls become causally connected 
as a whole and begin to contract, owing to their surface tension. This contraction ends up 
with the formation of a BH with a mass of the same order as that of the wall, when it enters 
the horizon. Therefore the size distribution of closed vacuum walls gets converted into the 
mass spectrum of resulting PBHs. Such a spectrum causes no problem for BBN predictions 
and satisfies all the other restrictions of possible PBHs contributions [12], while providing a 
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sufficient number of very massive BHs tliat could seed galactic nuclei^. The other generic feature 
of the mechanism [12] is the fractal structure of initial domain walls, which reproduces itself in 
a more and more decreasing scale. After the domain walls collapse, their fractal distribution, 
generated in fact during inflation gets imprinted into space distribution of less massive PBHs 
around the most massive ones. 

In the present paper we develop the ideas firstly revealed in [12, 13] and investigate the space 
correlation properties of the distribution of PBHs, which can modify the power spectrum and 
can lead to the fractal-like correlation properties observed in the galaxy distribution at least up 
to scales 100 Mpc [29]. The paper is organized in the following way. In the second section we 
describe the main issues of the mechanism [12] and discuss observational data in favor of early 
BH formation. The third section deals with the spatial distribution of the resulting clouds 
of PBHs. Finally, in section 4 we conclude and discuss possible effects of fractal-like PBH 
distribution. 

2 Formation of PBHs in inflationary PNG dynamics 

We review here a mechanism of primordial SMBHs [12] formation. This mechanism is based on 
the inflationary dynamics of a scalar fleld with two different vacuum states which get unequally 
populated after the end of inflation deeply in the Friedman-Robertson- Walker (FRW) epoch. 
Such a situation can occur quite naturally under the condition that the potential of the scalar 
flelds possess a valley along which the effective mass of the fleld vanishes comparable to the 
Hubble rate during inflation. The latter causes a large friction term in the equation of motion 
of the scalar fleld making it in fact to be effectively massless during inflation. The subsequent 
exit from the inflationary regime and further cooling of the Universe lead to decreasing of the 
expansion rate and therefore making the scalar field moving freely toward the minima of its 
potential. However it is the ability of quantum fiuctuations of effectively massless scalar field 
on the de-Sitter background [30] to redistribute the scalar field unequally between two different 
states. It makes the Universe to contain islands of less probable vacuum, surrounded by the 
sea of another, more preferable one. This means that once the vacuum states are populated 
at some time in FRW epoch, closed vacuum walls should appear in the Universe, surrounding 
the islands of less probable vacuum. When a closed domain wall becomes causally connected, 
it starts to collapse forming a BH with mass, being equal to that of the progeniting vacuum 
wall. It was shown in [12] that this mechanism can provide formation of sufficiently large closed 
vacuum walls, so that primordial SMBHs cold be formed with masses and abundance satisfying 
observations. 

Below we use the working example of a complex scalar field with the potential 

v{v)^ Kw? - fm' + moi (1) 

^Another possibility of massive PBH formation is considered in [27]; it involves the blue initial spectrum of 
fluctuation and subsequent accretion of quintessence. However, the last measurements of cosmological parame- 
ters by WMAP [28] disfavor the blue spectrum of initial perturbations. These measurements do not exclude the 
possibility of massive PBH formation by spikes on the background of the red spectrum, which can be realized 
in the mechanism [24]. 
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considered in [12], where 99 = r^e*^. This field is taken to be coexisting with the inflaton that 
drives the Hubble constant H during the infiationary stage. In this paper we accept the value 
H = 10^^ GeV for the Hubble constant during inflation, which is consistent with resent results 
of WMAP [25] on the energy scale of inflation. The term 

sv{e) = A'^{i-cose), (2) 

which reflects the contribution of instanton effects to the Lagrangian renormalization (see for 
example [31]), is negligible during the inflationary stage and during some period in the FRW 
expansion. In other words, the parameter A vanishes with respect to H. The omitted term 
(2) induces the effective mass of an angular fleld, and it begins to play a signiflcant role only 
at the moment, after inflation, when the Hubble parameter is equal to this effective mass in 
its gradual decrease with time {H — l/2t during the radiation-dominated epoch). Also, we 
assumed in [12] the mass rrir = \/A/ of the radial fleld component r^p to be always sufficiently 
large with respect to H, which means that the complex field is in the ground state even before 
the end of inflation. Since the term (2) is negligible during inflation, the fleld has the form 
if f/V^- e'^, the quantity f9 acquiring the meaning of a massless angular fleld. 

In [12] we applied the fact that the quantum fluctuations of the phase Sitter background 
can be pictured as a one-dimensional Brownian motion along the circular valley of the potential 
(1): the magnitude of 9 on the scale of the de Sitter horizon {H~^) fluctuates by an amount 

56 = H/27rf (3) 

per Hubble time {H~^). The amplitude of those fluctuations with wavelengths exceeding 
freezes out at its value because of the large friction term in the equation of motion of the 
phase, whereas its wavelength grows exponentially. Such a frozen fluctuation is equivalent to 
the appearance after every e-fold of a classical increment /decrement of the phase by factor 
66. Therefore assuming that the inflation begins in a single causally connected domain of size 
H~^, which contains some average value of phase 60 < tt one can conclude that more and more 
domains appear with time, in which the phase differs signiflcantly from the initial value ^o- A 
principally important point of [12] is the appearance of domains with phase 6 > tt. Appearing 
only after a certain period of time (e-foldings) during which the Universe exhibited exponential 
expansion, these domains turn out to be surrounded by a space with phase 6 < tt. 

The potential (1) possesses a U{1) symmetry, which is spontaneously broken, at least after 
some period of inflation. This spontaneous breakdown is provided by the condition of sufficiently 
large radial mass, = \l\>f > H. At the same time the condition 

m, = ^ «// (4) 

on the angular mass provides the freezing out of the phase distribution until some moment of the 
FRW epoch. After the violation of condition (4), the term (2) contributes signiflcantly to the 
potential (1) and explicitly breaks the continuous symmetry along the angular direction. Thus, 
potential (1) eventually has a number of discrete degenerate minima in the angular direction 
at the points 6min — 0, ±27r, ±47r, ... . As soon as the angular mass mo is of the order of 
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the Hubble rate, the phase starts be oscillating ^ about the potential minimum, initial values 
being different in various space domains. Moreover, in the domains where the initial phase is 
t: < 9 < 27r, the oscillations proceed around the potential minimum at Omm = 27r, whereas 
the phase in the surrounding space tends to a minimum at the point 9min — 0. This implies 
that closed vacuum walls Fig.l get formed ^ separating the island with 9min — 27r from those 
with 9„iin = 0. The surface energy density of a wall of width ~ 1/m ~ //A^ is of the order 
of ~ /A^. Those walls, which are larger than the cosmological horizon, still follow the general 
FRW expansion until the moment when they get causally connected as a whole; this happens 
as soon as the size of a wall becomes equal to the horizon size Rh ^. Evidently, internal stresses 
developed in the wall after crossing the horizon initiate processes tending to minimize the wall 
surface. This implies that the wall tends, first, to acquire a spherical shape and, second, to 
contract toward the centre ^. Should the wall at the same moment be localized within the 
gravitational radius, a PBH is formed. The phase distribution generated during infiation is 
reflected in the size distribution of vacuum walls and hence in the mass distribution of BHs 
appearing after the walls collapse. 

The above considerations are valid under the assumption that is situated at the vacuum 
value = f during the last 60 e-foldings of inflation, so that the only degree of freedom 
is the massless angular component f9. One can take into account the radial fluctuations 
which can actually kick r^^ out of the vacuum. The consequence of being kicked out of 
the vacuum over the top of the Mexican hat potential V{ip) is that strings get formed along 
the lines in space where r^p = 0. If it is the case the usual thermal scenario of the formation 
of the system of topological defects walls bounded by strings akin we have in the theory of 
invisible axion gets realized. Subsequently such structure of walls bounded by strings would 
disappear rapidly due to disintegration into separate fragments and futher emission of axion- 
like particles (see for review [32]). However the situation with radial fluctuations has been 
investigated both numerically and analytically, using the relevant Langeven and Fokker-Planck 
equations [26]. In particular it was found that r^p is kicked out of the vacuum if the Hubble 
constant during inflation exceeds the radius of the bottom of the Mexican hat potential H > f. 
In that case, strings form and the thermal scenario is reproduced. In some sense there is 
a "Hawking temperature" of order H during inflation which takes over to drive thermal-like 
phase transition. From the other side, as it was found in [26], sits in the vacuum whenever 
H < f during inflation. The last case is valid for the reasonable value of A ~ 10~^ and 
simplest exponential inflation. Therefore in the regime / lOH considered in present paper 
(see numerical parameters in Section 3) the only degree of freedom is the angular component 
f9, so only closed vacuum walls do form while the formation of strings is prevented. 

Generally speaking the domain walls we discuss, can decay by quantum nucleation of holes 
bounded by strings [36, 37]. Semiclassically, the hole nucleation is described by an instanton 

^Thc discussion on the coherent phase osciUations in the cosmology of the invisible axion see in [32] and 
references therein. Cosmology of rather heavy axion predicted in string theory is discussed in [33] . Variability 
of /, induced by inflaton, is discussed in [34]. 

•^The existence of such domain walls in the theory of invisible axion was first pointed out in [35]. 

''This condition can be violated for very large domains, as discussed in the next section. 

^The motion of closed vacuum walls has been derived analytically in [41, 42]. 
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and the decay probability can be expressed as [36, 37] 



P ^ A exp 



( 



~3^ 



) 



(5) 



where ^ ~ vr/^ and a ~ /A^ are the string and wall tensions, respectively, while the pre- 
exponent factor A can be calculated by analyzing small perturbations about the instanton [38]. 
Experessing (5) via parameters of (1) and (2) one can easily see that the probabihty to create 
a hole is suppressed by factor exp(— 10^(//A)^), which turns out to be very small owing to 
large hierarchy //A ~ 10^^ we use to realize the proposed mechanism (see the next section). 
Therefore on the timescale of the age of the Universe the walls are effectively absolutely stable 
relative to quantum nucleation of holes bounded by strings [36, 37]. 

There is also very unlikely that the initially closed walls can become open due to their mutual 
intersections. This argument can be confirmed by numerical analyzes [39], which demonstrated 
in the case of parallel sine-Gordon walls that such walls pass freely through one another. So 
it is unlikely that collisions of the walls leads to formation of holes and that closed walls 
can become open in such collisions. Due to non-perfect initial form of walls their flattening 
takes place in course of wall evolution. Such flattening changes the masses of the walls and, 
consequently, influences the mass distribution of black holes, which we discuss below. However, 
walls intersections in the course of flattening do not change the topology of walls, they do not 
convert closed walls into open ones bounded by strings. 

In the considered case the first (spontaneous symmetry breaking) phase transition takes 
place on the inflationary stage and neither angular nor radial fluctuations can provide formation 
of strings. It makes closed walls the only topological defect solution that can appear after 
manifest symmetry breaking, when the angular mass is switched on. Such walls are not bounded 
by strings by construction, and neither quantum hole nucleation nor mutual intersection can 
change the closed topology of these walls. Domain walls can become open when a string loop 
that binds a wall inside it hits a larger wall [37, 40]. Then pieces of two walls annihilate and 
a vacuum hole bounded by string is produced. However as we argued above all the walls are 
initially closed and have no strings on their boundary, so during intersections they will pass 
each other without annihilation [39] thus preserving their closed character. 



The general conditions for the formation of primordial black holes have been presented in [13]. 
In that paper we also performed calculations for formation of walls of smaller sizes inside the 
wall of larger size. The conclusion is that their total mass is only twice a mass of main wall, 
which determines the mass of the largest black hole. BHs are formed as the result of the closed 
walls evolution. Hence, total mass of black holes disposed closely to the largest BH increases 
weakly its effective mass. 

The present section contains calculation of BH spectrum around the main BH. The size 
of area which we consider is much greater (in orders of magnitude) than the size of the wall 
progeniting the largest BH. It contains a lot of smaller walls what strongly increase total mass 



3 Space distribution of PBHs 
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of the BH cluster. So that is of great interest to study the properties of distribution of BHs 
with intermediate masses around the most massive ones over those scales which can be roughly 
attributed to the galaxy clusters scale. 

General strategy is as follows: we choose the scale of some region containing BH with the 
given value of mass; then we consider what the value of phase 6^ must have been when this 
region was in a single horizon. We do that by deciding what value of 9 would have (on average) 
led to a single black hole of the given mass. Then we use that value of Or to find the density of 
BHs of smaller sizes in the considered region. 

Let Mo be the mass of most massive BH and M the integrated mass of the coexisting BHs 
inside a sphere of radius r centered at the most massive BH. It is convenient in our mechanism 
to trace the space distribution of resulting BHs by the distribution in size of their progenitive 
closed domain walls. In this sense, there are three scales related to our problem, namely the 
size of the closed wall i?o, which finally gets converted into the most massive BH; a size R of 
a closed wall that becomes a BH of intermediate mass inside the above-mentioned sphere; and 
the radius r of the sphere itself. Every mentioned scale can be attributed to the number of 
e-foldings left until the end of inflation. We assign these numbers of e-foldings as A'^q, and 
Nr respectively. It is clear that a region of size r formed at the inflationary stage will follow 
the expansion of the Universe, depending upon the current equation of state, and hence we can 
approximately admit 

r = r(7V„0=Mt)e^'-^-, (6) 

where ruit) is the size of the Universe at time t. Let us assume that the angular component 
of our field has the phase 6^1 once a region of size r leaves the horizon when the Universe still 
has Nr e-foldings to inflate. So, at some later time, characterized by Nq < Nr e-foldings before 
the end of inflation, one has already e^'^^''^^''^ causally disconnected regions inside the volume 
created at e-folding Nr. As we saw in the previous section, the phase 9 fluctuates at each 
e-folding like a one-dimensional Brownian motion. According to our mechanism, the biggest 
domain wall is formed once the phase 9 overcomes tt for the first time, being initially (at the 
beginning of inflation) somewhere in the range [0; vr] (see Fig. 1). It is this wall that collapses 
into the most massive BH serving as the centre of BHs cluster. By assumption, only one of 
the above-mentioned domains contains a phase larger than tt at e-folding number Nq. So the 
probability of finding the phase within the interval [tt; 7r±(5^] in one of the causally disconnected 
regions at e-folding number A^o is equal to l/e^^^''~-^°\ 
On the other hand, the well known formula [43, 44, 45] 



W(N, 9) = , exp 
^271 {Nr - N) 



(Or - 9 f 



(7) 



represents the probability to find phase 9 in the interval [9; 9 ±69] at arbitrary e-folding number 
within causally disconnected volumes with the size of order oil/ H provided they had unique 
phase 9r at e-folding number A^^- The value 59 is given by formula (3). Thus the number of 
causally independent domains acquiring the phase in the interval [tt; tt ± 5^] just during the 
e-folding number N is 
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X(Ar) = W^(Ar,7r)e3(^'-^). (8) 

The value of the phase Of. can be expressed by equating number of domains containing the 
phase in the interval [tt; tt ± 50] at e-folding number Nq to unity (by definition) 

X(iVo) = 1, 

what gives 



Or^-K- Se^{Nr - No) [6 {Nr - Nq) - ln(27r {Nr - No))]. (9) 

Combining Eqs. (7), (8), (9) the probability of formation of walls with smaller masses at 
e-folding number N can be obtained as follows: 



[2% {Nr - No 



2(Nr-N) 

K (N) = ^'"'-^ — exp 

/27r {Nr - No) 



3{Nr-No) 
Nr-N 



2 



+ 3{Nr - N) 



(10) 



The formulae written above are formally valid for arbitrary parameter A'^o which is connected 
with the mass of main BH. The BH of any mass produces cluster of smaller BHs in the manner 
discussed above. The choice of the mass of main BH could be done using mass distribution 
of BHs calculated earlier - see Fig. 2. One can see from the figure, that amount of BHs with 
mass 5 x lO^M© - see Fig. 2- is about 10^^, i.e. approximately coincides with the total number 
of galaxies. To treat the BH cluster structure as the pattern for future galaxy formation we 
choose below the parameter A'^q, corresponding to this value of mass for the main BH. 

Now we shortly remind the connection between mass of black holes and e-folding accord- 
ing to [12]. The mass M of a closed wall is connected with its size 

M « 47ri? (AT) V, (11) 

where a is the surface energy density of the wall. 

There are two different scales related to the spherical wall dynamics. The first one Ri{N) 
is determined from the condition that the wall mass (11) is equal to the energy of the plasma 
contained in the domain bounded by this closed wall at t = R, Ey = p^R^, where p = ^g*T^ 
is the energy density of the plasma with g* degrees of freedom during the radiation-dominated 
epoch. Taking into account that 

a = 4AV, (12) 
and equating energies Ey and M from expression (11), one obtains the scale 

R-j, (13) 

which also can be expressed through the temperature T when the above-mentioned equality 
was established 

R = Lu^e^-^^; (14) 
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here Tq and Lu are temperature and size of the modern Universe. Combining the formulas 
written above, one obtains the maximal size of a spherical wall at which it does not influence 
the law of local expansion, when it enters horizon: 

Rm - i^^g) {LuTof el(^-^-). (15) 



2.7 X 10 



29 



Another scale arises if we take into account the internal pressure of the wall. Its influence 
is developed after the wall comes under the horizon and leads to the termination of the wall 
growing and its acquiring a spherical shape. This moment is characterized by the instant t2 
such that 

lhorit2) = R{t2) = i?2, (16) 

where the size of the horizon during the radiation-dominated era scales like lhor{t2) = 2^2- We 
also recall that walls get nucleated under the violation of condition (4), which happens once 
the friction term ?>H9 in the equation of motion for the phase 9 becomes comparable with the 
'force' term m^^, i.e. at the moment t* when 

H{t*)^me. (17) 

Taking into account the expression for the Hubble parameter H — ^J^j^, and the connection 

of the temperature with time T = ( ^^^1^ ) ^ ^/^, one obtains the time moment and the 

1/4 



temperature of walls nucleation in the early Universe 



Applying these formulas to (16) and taking into account that the seeking scale can be expressed 
in the same manner as (14), with T* instead of T, one can easily flnd the size of wall when it 
enters horizon 

MN)^\{^)\-'^-^^K C^(j^)""^. (19) 

Thus the maximal comoving size of a spherical wall at FRW epoch seeded at e-folding number 
N, is determined from the condition 

R{N) = min(i?i(iV), R2{N)). (20) 

The condition (20) constrains the size of a spherical wall so that its contribution into the local 
density does not dominate when it enters cosmological horizon. 

Indeed, for a sufficiently big wall, of radius exceeding Riyj^ = given by the constraint 

(15), its pieces start to dominate over radiation energy within the cosmological horizon, while 
the wall is still causally disconnected as a whole. Therefore the local expansion rate in a causally 
connected area dominated by a piece of the wall turns to be aioc{t) oc t^. Such superluminal 
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expansion makes the successive evolution of this region elusive for the distant observer. It 
puts the upper limit (15) on the maximal mass of BHs that can be formed by the considered 
mechanism from a spherically symmetric wall. However account of the non-perfect initial form 
of the wall and of the effect of crinkle flattening, considered in the next section, offers the 
possibihty to avoid the constraint (15) even for supermassive walls. Then the condition (20) 
can be relaxed and the maximal size of wall can always be given by the size it has, when enters 
the horizon. Local crinkle decay dominance (see the next section) leads, however, to deviation 
of local expansion law from the pure RD case, what modifies the definition (19) for the maximal 
size of wall. 

According to the mechanism under consideration, the energy of the wall with surrounding 
particles E^, will be converted into a BH of mass: 

Mbh{N) « (21) 

This happens when the gravitational radius Vgrav — '^Eyj/Mp is larger than the minimal size 
of contraction of the collapsing wall. This minimal size of contraction does not exceed a few 
times the wall's width. 

Using formulas (6), (8) and (21) one can calculate the desired distribution of intermediate 
BHs in the vicinity of the most massive one. 

In Fig. 2 following [13] we present the mass spectrum of BHs produced by the described 
mechanism. Such spectrum strongly depends on the parameters of the underlying field model. 
Here these parameters are chosen in such a way that the average number of BHs with mass 
~ 10^ Mq, is about 10^^ what coincides approximately with a number of galaxies in the visible 
Universe. It makes the difference of the model in Fig. 2 from the results of [13]. 

The mass distribution of intermediate BHs around the most massive one is presented in Fig. 

3 for the same parameters of the potentials (1) and (2) as have been chosen to produce the 
distribution Fig. 2. Actually we always deal with a cluster of BHs centered by the most massive 
one. This BHs cluster contains from the beginning ~ 0.3 x 10^^ BHs of smallest masses about 
~ IO^^Mq and massive BH with mass ~ 2.7 x 10^ Mq in its center. The total mass of the BHs 
in the cluster is ~ 0.3 x lO^^M© and is comparable with the mass of a baryon component in a 
galaxy. Analysis of Fig. 3 reveals that average density decreases quickly with size growing. For 
example, average density at the scale 1000 parsec in 10^ times smaller comparing with average 
density in the center, at the scale of order lOpc. 

4 Possible role of crinkles 

Let us consider the problem of squaring of a wall surface just after the wall formation. By 
definition, the spatial position of any wall is determined by the surface with phase value 9 = tt. 
Let we have an approximately fiat surface determined as space points containing the phase 
^ = TT of size 1 / ifg at the e-folding number A^. Then, one e-folding later, the surface will be 
divided into 1/e^ parts. Approximately half of them embrace new value of phase 9 > ir while 
the remaining half acquires the values of phase 9 < n. The total surface of interest, with the 
phase 9 — TT, appears to be warped and increases by a factor q varying in the range ~ 1-2. The 
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following cosmological expansion increases the sizes of crinkles exactly like total size of a closed 
wall. If this process lasts for, say, N e-foldings, the surface increases by a factor in addition 
to the standard inflationary expansion by a factor . 

Here some explanation is necessary: In the course of inflation the warping goes to smaller 
scales, so that smoothing of large crinkles beyond the horizon is accompanied by creation 
of small crinkles within it. It should be reminded that in this process not walls, but physical 
conditions for their origin are created. This fractal-like process of crinkles formation terminates 
when the size of the wall becomes comparable with the wall width d = 1/mg ^ 0. The crinkles 
of sizes smaller than the wall width d will be smoothed out during the period of wall formation 
and their effect could be ignored in the analysis of successive wall evolution. Hence our direct 
aim is to determine the e-folding at which crinkles of size d are still produced. 

The scale of the crinkles produced at e-folding number is approximately equal to the 
wall width d 

d ^ Icrrnkle = rc;(f )e^<^-^- = Lc;^e^<^-^-. (22) 

This formula immediately leads to the number of e-foldings we need: 

(23) 

Therefore the mass of a wall just after its nucleation is given by the expression 

^ aAnR (Nf c;^-'^^. (24) 

Choosing the realistic numbers N — = 10-20 and ^ = 1.5, one can see that the effective 
mass of the wall increases at the moment of its formation by a factor ~ 10^-10^. Also one 
has to note that the size Ri{N) we have calculated in the previous section should be initially 
corrected by a factor once the effects of crinkles are taken into account. 

The internal tension in the wall leads to its local flattening with the growth of the horizon. 
The released energy is converted into the kinetic energy of the wall and into an energy of 
outgoing waves, which are predominantly the waves of that scalar fleld that produces the 
domain walls. In our case, it is the phase 9 that acquires dynamical sense of the angular 
component of fleld ip. 

A typical momentum k of the waves is of the order of the inverse wall width 1/d, hence 
k ~ nig and the waves are necessary semirelativistic. Their kinetic energy decreases due to 
the cosmological expansion. So particles of mass mg are concentrated in the vicinity of the 
wall. These particles can further contribute to the process of BH formation. The situation can 
become more complex if the particles are unstable. In this case, the products of their decay 
are relativistic and most of the time, they escape capture by the gravitational fleld of the wall. 
The main immediate effect of these relativistic products is the local heating of the medium, 
surrounding the region of BH formation. Their concentration in this region and successive 
derelativization can later provide their possible contribution to the formation of BHs and BH 
clusters. 



Na^Nu + In 



T*d 
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The density of these particles is diluted by cosmological expansion, whereas the wall kinetic 
energy is mainly redshifted, making the expression (24) approximately valid with ~ 1 in the 
period of BH formation, when the whole wall enters the horizon. 

However, the contribution of crinkles can be significant, in the period when the wall enters 
the horizon, and its successive evolution proceeds so quickly that the effects of wall flattening 
are not suppressed by the cosmological expansion. In this period the wall flattening and crinkle 
decays transform the energy, stored in the smaller-scale crinkles, into the energy of scalar field 
waves and the kinetic energy of the wall. These products of wall flattening dominate in the 
local energy density, thus preventing the superluminal expansion of the considered region. If 
the effect of crinkles is sufficiently strong, the wall dominance may not arise, since the energy 
released in the fiattening is always dominant over that of the negative pressure of the fiattcned 
wall within the horizon. The RD (or MD) regime of dominance of the energy, released in the 
fiattening, implies the difference between the expansion law within the considered region and 
the general expansion; however, being subluminal, it continues until the whole wall enters the 
horizon. 

The longer the period of crinkle-decay-products dominance is, the larger is the overdensity 
in this region, as compared with the mean cosmological density, (Sp/p) ^ 1, when the wall 
enters the horizon. This difference provides the separation of this region from the general 
expansion and the effective BH formation in it, similar by to the PBH formation on the RD 
stage, considered by [18] (for a review see [17, 23]). This effects can modify the numerical values 
in BH distributions presented in the previous section. A special study is, however, necessary 
to take it into account. 

5 Discussion 

It should be obvious that the mechanism proposed in this article, which explains the formation 
of PBHs by the pregeneration of objects that build up the mass of SMBHs in the centre of the 
galaxy, works in reality by production of BHs with smaller masses that accompany the most 
massive BH. An indirect hint [47] that a situation of this sort can arise has been found recently 
in the centre a rather famous blazar, Mkr 501. Namely the combination of observations of the 
host galaxy allowing a direct measurement of the mass of the central BH in Mkr 501 and those 
of the modulation of high energy emission might be brought into agreement by assuming a 
binary system of highly unequal SMBHs hosted in the centre of this blazar. Thus, having BHs 
of intermediate masses accompanying the most massive ones, which are capable of coupling 
into binary systems, turns out not to be unrealistic. 

We investigated in detail the primordial spatial distribution of BHs with masses in the 
range ~ 10~^-10^Mq. It was shown that such BHs could concentrate around a massive one 
in the scale of up to 10 Mpc. The total mass of a group of these primordial BHs can be of 
the order of a mass of ordinary galaxy clusters. The number of such a groups could be made 
to coincide approximately with a number of galaxy clusters. Their formation begins at a very 
early stage, which is characterized by the temperature T* ~ yjMpme ~ lO^GeV; for the most 
massive BHs, it may continue even after the period of Big Bang nucleosynthesis, but stop in 
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any case long before star formation. Being compatible with the constraints on the average 
primordial light element abundance and with the averaged restrictions on the distortions of 
the CMB black-body spectrum, the formation of PBH clusters around the most massive BH 
may lead to local peculiarities of pregalactic chemical composition, and to specific effects in 
the spectrum, angular distribution and polarization of small-scale CMB fiuctuations, accessible 
to the observational test. After recombination of hydrogen at the modern matter-dominated 
stage, baryons are influenced by the gravitational potential of a primordial BH group, which 
facilitates galaxy formation. The proposed model can offer the alternative to the standard 
picture of galaxy formation by the slow growth of baryon-density fiuctuations. 

To illustrate our approach we considered the simple model of a PNG field in an infiationary 
Universe. The physically self-consistent realization of the proposed scenario should imply more 
sophisticated models. Such a realization should provide the compatibility of this scenario with 
the observational constraints. Namely, the large-scale fluctuations of the PNG fleld should 
be suppressed so as not to cause too large CMB anisotropy. Another problem of a realistic 
scenario will be to satisfy the constraint on the oscillations of the energy density of the PNG 
field, playing the role of cold dark matter, or to provide a sufficiently rapid decay of PNG 
particles, if the parameters of the chosen model lead to the contradiction with this constraint. 

At the end of the discussion it should be noted that the existence of such systems of pri- 
mordial BHs correlated in space offers a new element to the standard lambda-CDM model. It 
could add additional power to the spectrum of density fiuctuations on small scales, where the 
measurements of the Lya forest are relevant [48] . It is important that the described structures 
of PBHs do not evaporate significantly due to the large amount of their constituents. This fact 
makes the clusters of PBHs very favorable to being gravitationally bounded during the epoch 
of galaxy formation, up to the galaxy clusters scales. 
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Figure 1: The explicit deconvolution of the potential (2) (top) shows that the vacuum manifold 
of angular component 6 consists of two stable states associated with and 27i akin to the simple 
double-well potential (only the two closest minima are shown). This situation is an analogue of 
that with the breaking of discrete symmetry. A domain wall occurs at the boundaries, between 
regions of space with the vacuum value of the phase 9 in different states, with 9 interpolating 
between and 27r across the wall. The energy of the configuration is rising over the potential hill. 
Two neighboring domains (bottom) filled with the phase values above and below tt, corresponding 
to the two different vacuum states, give rise to the formation of a domain wall once the motion 
of the angular component becomes underdumped, driving the phase to different minima, namely 
27r and respectively. 
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Figure 2: Explicit mass (in g.) distribution of PBHs, simulated for the following specific pa- 
rameter: A = 100 GeV, / = 8 X 10^^ GeV. The total mass of BHs over the Universe is IQ^^g. 




Figure 3: Total mass of black holes within the sphere of radius r (in pc). Parameters of the 
model: A = 100 GeF, / = 8 x 10^^ GeV; the mass of the central BH is 2.7 x lO^M© 
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